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Implications of the unusual structure in the pp correlation from Pb+Pb collisions at 

158 AGeV 
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The recent NA49 measurement of two-proton correlation function shows an interesting and un- 
expected structure at large relative momentum. Applying source imaging techniques to the mea- 
surement, we find an unusually steep drop-off in the two-proton source function. We show that 
the steep drop-off is due to the structure in the correlation and the drop-off cannot be explained 
using conventional correlation analysis. We suggest possible physics reasons for the unusual source 
function. 

PACS number(s): 25.75.-q, 25.75.Gz 

The NA49 collaboration has recently measured the 
two-proton correlation function from Pb-|-Pb collisions 
at 158 AGeV, integrated over transverse momentum in 
the rapidity range 2A < y < 3.4 (midrapidity « 2.9) [Q. 
The measurement shows an interesting structure around 
q'inv (the magnitude of the proton momentum in the pair 
center-of-mass frame) of 70 MeV/c, which statistically 
significantly deviates from the expected unity. NA49 has 
extensively studied possible systematic effects and none 
of these systematics can account for the structure. In 
this letter, we pursue possible physics explanations for 
the structure. 

Similar structure has not been seen in other experi- 
ments for different collision systems or energies. Sim- 
ilar structures have been predicted in the tails of two- 
pion correlations and have been ascribed to various 
effects, from hidden correlations to a breakdown of 
the smoothness approximation Applying these pre- 
dictions to two-proton correlations is difficult because of 
the differences between the two-proton and two-pion final 
state interactions. Source imaging is an ideal tool 
to separate effects due to two-proton final state inter- 
actions and wavefunction antisymmetrization from those 
due to the source function itself. Moreover, one can ob- 
tain the source function in a model-independent manner, 
i.e. without a Gaussian source assumption. 

We have inverted the NA49 two-proton correlation 
function and obtained a two-proton source function that 
is exceptional in several ways. Our naive expectation was 
that the source would appear somewhat Gaussian (as is 
the case for all two-pion and two-proton sources that have 
been obtained so far @-|^). Instead, the source is consis- 
tent with a step-function with radius 10 fm. In some 



sense this is not a surprise as the structure in the tail 
of the NA49 correlation function is reminiscent of the 
Fourier transform of a sharp object. In fact, the edge 
can be removed by smoothing out the structure. As if 
finding a step-function source is not strange enough, this 
source cannot be explained within the standard Koonin- 
Pratt formalism: in the Koonin-Pratt formalism, the 
Fourier transform of a source must be positive every- 
where while the Fourier transform of the NA49 source is 
not. Thus, using imaging and this test of positivity, one 
may infer whether a given source is breaking something. 
Knowing this, the task then becomes understanding what 
is broken. 

The outline of this letter is as follows. First, we de- 
scribe the conventional expectations for sources in the 
Koonin-Pratt formalism. Second, we detail why the 
Fourier transform of a source in the Koonin-Pratt for- 
malism must be positive, thus giving us a test to see if 
a given source fits in this formalism. Then, we perform 
the inversion of the NA49 two-proton correlation function 
and show that it fails this test. We demonstrate that the 
structure in the tail of the correlation function is the ori- 
gin of the failure. Finally, we discuss physics that might 
be missing from the Koonin-Pratt formalism that could 
account for the unusual behavior of the source. 

We begin with our expectations for the source function 
in the Koonin-Pratt formalism |§-|l^. A two-particle 
correlation function can be written as a convolution of a 
source function, Sq' (r ), with the Wigner transform of the 
particle pair wavefunction, g^{r, q'), in the pair center-of- 
mass system (c.m.s.) Q: 



C{q) =Jd^ 



d^q' 
(27r)4 



9qir,q)Sg>{r). 
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Here q = ^{pi ~ P2) is the asymptotic particle momen- 
tum, q' is the internal four-momentum (and need not cor- 
respond to on-shell particle momenta), and r = fi — r2 
is the separation of emission points. Both (7g(r, q') and 
q' are defined through 



g^{T, q') = j e<-^'^{r + C/2)<I>,-( 



r-C/2). (2) 



Here ^q{r) = (/)g'(r )e^**'"" is the particle pair wavefunc- 
tion. Note that the different phases of the wavefunction 
in Eq. (|^) make independent of the relative emission 
time. 

In the Koonin-Pratt formalism, one makes two as- 
sumptions. The first is that the source function only has 
a weak dependence on relative momentum correspond- 
ing to off-shell particles, i.e. Sq'{r) « S^{f). This is 
a smoothness assumption [p|Jl2|. With this, the integral 
over the internal momentum q' in Eq. (^) can be per- 
formed and we obtain 



C((f)= / d'r\Mr)\ Sqif). 



(3) 



One further assumes that the two-particle source is a con- 
volution of single-particle sources (i.e. particles are emit- 
ted independently): 

%(r ) = 

(4) 

where D is the normalized single-particle source in the 
c.m.s. and has the conventional interpretation as the nor- 
malized production rate of final-state particles. The inte- 
gral of D over the particle emission times, / dtD{r, t; q), 
is then the normalized phase-space distribution of parti- 
cles after last collision (freeze-out) and S^{f) is the prob- 
ability density for a pair each with momentum q to be 
emitted a distance f apart in the c.m.s. 

In Eq. (^) we do not need to consider the contribution 
to the source from pairs with large relative momentum 
{\q \ = q > gcut « 100 MeV/c): the kernel cuts off the 
contribution from these pairs. This happens because the 
source varies slowly on the length scale of the oscillations 
in the kernel for q ^ gcut, hence the integral in Eq. (|^) av- 
erages to zero. The fact that only low-g pairs contribute 
to the correlation is used to justify dropping the q depen- 
dence from the single-particle sources entirely ]lO|Jl^ . 

The Fourier transform of Eq. (0) is 



S{k)= / d^re'^-'^Sgif) = / d^Rd\e'^''x 



dtiD[R+^,ti;q] dt2D(R-^,t2;~q 



b{k,q) D*{k,-q), 



where D is the Fourier transform of the time-integrated 
single-particle source. // we can neglect the q dependence 
of the source function (or equivalently the momentum 
dependence of the single-particle sources in the c.m.s.), 
then we find the following interesting result: 



~S{k) 



d\e'^-''S{f) = \D{k)f > 0. 



(6) 



So, checking the Fourier transform of the source for posi- 
tivity can serve as a test of the validity of the underlying 
assumptions. 

For pions, we note that Eq. (||) is actually a restate- 
ment of the theorem that the correlation function must 
always be larger than one (after the Coulomb correc- 
tion). To see this, first realize that sources of indis- 
tinguishable particles must be symmetric under coordi- 
nate reversion: r — > —r. Using this, Eq. (|^) becomes 
S(k) = J d^r cos(r • k)S{r) > 0. Now, for like-pions 
Eq. (|) gives C{q) - 1 = J d^r cos(2f • q)S{f) = S{2q) 
which must be positive. 

We now turn to the practical issues of extracting in- 
formation about the source from the NA49 correlation 
function The NA49 correlation function is averaged 
over the orientation of q so we integrate out the angular 
dependence of the correlation, arriving at 



C{q) -1 = At: j dr r"^ K{q,r)S^{r). 



(7) 



Here q = \q\ and r — \r\. The angle averaged kernel, 
K{q,r) — J d{cos6^f) \(j)^{r )\'^ — 1, can be written as a 
sum over the two-proton partial waves: 



K{q,r)^ I + 



9%{r) 



1. 



(8) 



In what follows, we calculate the two-proton relative 
wavefunctions by solving the Schrodinger equation with 
the REID93 and Coulomb potentials for partial 
waves with £ < 2 and use pure Coulomb waves for those 
with £ > 3. Note that the results are not noticeably al- 
tered by using the full solutions to the Schrodinger equa- 
tion for the higher partial waves due to the low relative 
momentum of the pair. 

Since both Eqs. (^) and (0) are simple integral equa- 
tions with non-singular kernels, they may be inverted 
to obtain the source function. We begin by discretiz- 
ing either Eq. (|^) or (Q) to obtain the matrix equation 
Ci — Kij Sj . If the data were measured with infi- 
nite precision we could invert this matrix equation to 
find Sj = Since we must account for the 

finite measurement errors, we proceed as in j^-Q] and 
find the set of source points, Sj, that minimize the x^. 
Here, = - K.jS,)^ A^C,. This source is 

Sj = Y.i{iK^BK)-^K^B]ji{C^ - 1) where is the 
transpose of the kernel matrix and B is the inverse co- 
(5) variance matrix of the data -By =Sij/A'^Ci. The error on 
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the source is the square-root of the diagonal elements of 
the covariance matrix of the source, A^S* = {K^ BK)~^ . 

In this inversion process, all q dependence of the cor- 
relation is ascribed to the kernel. If only the kernel has 
a q dependence, then the imaging of Eq. (^) or (0) is 
unique. On the other hand, if the source also has a q 
dependence, then the imaging gives the source averaged 
over q. While only pairs with q ^ gcut contribute to the 
correlation, this does raise the possibility that a strong q 
dependence might cause trouble with the imaging. This 
possibility will be investigated in a future article [|l5|. 

Averaging both sides of Eq. (^) over angle, we obtain a 
condition for the angle averaged source (which we assume 
to have no q dependence): 
1*00 

S{k) =4:71 dr r^S{r)jo{kr) > 0, (9) 
Jq 

where jo{kr) is the 0*'' order spherical Bessel function. 
The Fourier transform is a linear operation, so the errors 
propagate via AS{k) = in dr r^AS{r)jo{kr). 

We now turn to the NA49 two-proton correlation 
data In Fig. |^a, we reproduce the data with the 
filled circles. Note the oscillation at q ~70 MeV/c and 
maybe at ^^100 MeV/c. These oscillations are respon- 
sible for all the interesting behavior we discuss below. 
Figure |^a shows the source function obtained by imag- 
ing the correlation function in Fig. No constraints 
have been imposed in the imaging. Note that the source 
is nearly flat up to r ~ 10 fm, after which the source 
drops rapidly to zero in the bin centered at 14 fm. The 
source remains several orders of magnitude smaller than 
the main part of the source out to the edge of the image 
at r w 65 fm. For comparison, we plot NA49's best-fit 
source in the dashed curve in Fig. 2a. NA49 fit a Gaus- 
sian single-particle source to the q < 48 MeV/c points 
of the correlation function and found a best-fit Gaus- 
sian radius of 3.85 fm. Examining Fig. 2a., the best-fit 
Gaussian source describes the r < 10 fm points and fails 
beyond that. As a side benefit of imaging, we can ex- 
tract the integral of the source over the imaged region. 
We find 47r J^^^"" dr r^S{r) = 0.78 ± 0.55. We should 
mention that the constant weak decay correction applied 
in ref. § does increase the magnitude, and hence the 
integral, of the source over the uncorrected one, but the 
correction does not change the structure of the oscilla- 
tion. 

As a consistency check, we uninverted the source into 
the correlation function which is shown as the solid his- 
togram in Fig. The resulting correlation function is 
consistent with the original but the structure in the tail 
of the correlation has been washed out. This is easily un- 
derstood as an effect of the loss of information in going 
from a correlation function with ~ 20 independent data 
points to a source with only 7 points. 

Now we test whether the imaged source fails the 
Fourier transform test discussed above. We must take 



care in doing the integral in Eq. as simply discretiz- 
ing it introduces numerical error at high fc, exactly where 
a failure of the Fourier transform test will occur. Instead, 
we use an interpolation procedure from ref. [l^ to eval- 
uate Eq. (|) with the S{r) in Fig. §i. The result of this 
integration is shown in Fig. |^a. S{k) goes negative in the 
range 100-150 MeV/c. 

For comparison, in Fig. |^b. we replace the structure 
in the tail of the correlation with statistical noise. To 
be exact, we replace the measured correlation function 
in 36 < (7 < 102 MeV/c with a correlation consistent 
with one, but containing the original errors and statis- 
tical scatter. Applying the imaging technique, we ob- 
tain the source function shown in Fig. ^p. For compari- 
son, we plot the best-fit Gaussian source along with the 
"smoothed" source in Fig. 2b and the two sources are 
consistent. Furthermore, not only does this modified 
source appear less flat and somewhat Gaussian but its 
Fourier transform is positive as shown in Fig. |^b. Thus 
we have demonstrated that the structure in the tail of the 
correlation is the cause of the failure shown in Fig. ||a. 
This assertion should not come as a surprise as the be- 
havior of the correlation at the scale q is responsible for 
source structure on the scale r ^ T:hc/q. 

At this point we comment on the model comparisons 
by NA49 using different single-particle and two-particle 
sources [|l|. First, NA49 compared their data to corre- 
lations generated from VENUS and RQMD simulations. 
Both models generate the final freeze-out distribution of 
protons (equivalent to our _D(r, i;g)) and these distribu- 
tions are input to a correlation afterburner that performs 
the integrals in Eqs. (^) and Thus, by construction 
these correlations pass the Fourier transform test and 
cannot describe the tail of the correlation data. NA49 
also tried a hard edged two-proton source but could not 
reproduce the observed structure in the correlation with- 
out abandoning the independent source assumption. 

We are now left with a dilemma: something is wrong 
with either the imaging, the assumptions in the Koonin- 
Pratt formalism or both. While experimental confirma- 
tion of this new data is needed, we suggest three possible 
physics sources of the problem (as we have indicated in 
italic text in this letter): (1) a strong low-q dependence 
of the source that is not imaged, (2) the source not being 
a convolution of single-particle sources, or (3) a break- 
down of the smoothness assumption. 

In the case that there is a strong low-g dependence 
in the NA49 source function, such a q dependence can- 
not be reconstructed by the imaging. In this case, the 
Fourier transform test is inappropriate as the imaged 
source does not refiect the true NA49 source. One would 
think that a strong low-q dependence of the source might 
result from strong position-momentum correlations in the 
single-particle sources. In ref. [|l^, the authors delib- 
erately increased the degree of position-momentum cor- 
relation by aligning the proton transverse momentum 
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with the transverse position of the proton at freeze- 
out in RQMD. Even after this dramatic increase in the 
position-momentum correlation, they were able to image 
the source function reliably. 

In the second case, the NA49 two-proton source is not 
a convolution of two independent sources. In this case, 
both imaging and the Fourier test are appropriate. How- 
ever, interpreting the source function becomes difhcult 
as we cannot model it simply with a transport model. 
Makhlin and Surdutovich Q suggest that this possibility 
could occur if there are hidden correlations in the system. 
Such hidden correlations would arise from having three 
or more dynamically or statistically correlated particles 
in the final state while only observing two of them. 

Finally, in the event that the smoothness assump- 
tion is no longer valid, the entire formalism breaks 
down as we can no longer perform the integral over the 
internal four- momentum in Eq. (|^). This is the worst 
possible case from our standpoint as neither Eq. ^) nor 
Eq. (0) is valid and we cannot apply our imaging pro- 
cedure. Pratt 1^ argues that this could be caused by 
strong off-shell effects and these effects would be mag- 
nified by strong position-momentum correlations in the 
source and by the short-range interaction between the 
protons. Pratt comments that these effects are larger 
in small systems, e.g. with single-particle source size 
R 1 fm, much smaller than the source imaged here. 

In conclusion, we have imaged the recently measured 
NA49 two-proton correlation function and found an ex- 
traordinary two-proton source function. This source 
function appears to be a step function and we show that 
the steep fall-off in the source is due to the structure in 
the tail of the NA49 correlation function. Furthermore, 
we show that this source does not fit into the Koonin- 
Pratt formalism. While we have suggested several pos- 
sible physics explanations for the structure in the two- 
proton correlation, its origin remains a puzzle. 
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FIG. 1. (a) Correlation function from ref. jl] (filled circles, small caps on error bars) and the restored correlation function 
(solid histogram, large caps on error bars), (b) Same as (a) but the structure between 36 < 5 < 102 MeV/c is replaced by a 
fiat correlation with experimental error and realistic statistical scatter (filled circles, small caps on error bars), and the restored 
correlation function (solid histogram, large caps on error bars). 
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FIG. 2. (a) Source function obtained by imaging the correlation data (solid histogram) and NA49's best-fit Gaussian source 
(dashed curve), (b) Source function obtained by imaging the "smoothed" correlation function (solid histogram) also with 
NA49's best-fit Gaussian source (dashed curve). 
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FIG. 3. (a) Fourier transform of the source function in Fig. ga. Notice the dip below zero in the range 100 < k < 150 MeV/c. 
(b) Fourier transform of the source function in Fig. Notice that this function is everywhere positive. In both panels, the 
vertical lines represent the error band of the calculation. 
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